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We present a time-lapse sequence of 20 near-infrared (J- and X'-band) snap- 
shots of the central 20" x 20" of the Crab pulsar wind nebula, taken at sub- 
arcsecond resolution with the Hokupa'a/QUIRC adaptive optics camera on the 
Gemini North Telescope, and sampled at intervals of 10 minutes and 24 hours. 
It is observed that the equatorial wisps and polar knots in the termination shock 
of the pulsar wind appear to fluctuate in brightness on kilosecond time-scales. 
Maximum flux variations of ±24 ± 4 and ±14 ± 4 per cent relative to the mean 
(in 1.2 ks) are measured for the wisps and knots respectively, with greatest statis- 
tical significance in J band where the nebula background is less prominent. The 
J and K' flux densities imply different near-infrared spectra for the nonthermal 
continuum emission from the wisps and outermost polar knot ('sprite'), giving 
F v oc //-°- 56±012 a nd F u oc //-°- 21 ± 013 respectively. The data are compared with 
existing optical and UV photometry and applied to constrain theories of the vari- 
ability of the wisps (relativistic ion-cyclotron instability) and knots (relativistic 
fire hose instability). 

Subject headings: ISM: individual (Crab Nebula) — ISM: jets and outflows - 
pulsars: individual (Crab Pulsar) — stars: neutron — supernova remnants 



1. Introduction 

Many young pulsars in supernova remnants are embedded in synchrotron nebulae, 
known as pulsar wind nebulae (PWN), containing relativistic electrons and magnetic flux 
emitted by the central object. Multiwavelength imaging at subarcsecond resolution reveals 
that many PWN share a common morphology, consisting of (i) bipolar jets of unequal bright- 
ness, directed along the (inferred) pulsar spin axis and terminated by one or more bright 
knots, and (ii) fibrous arcs, or 'wisps', that are cylindrically symmetric about the spin axis, 
concave or convex with respect to the pulsar, and associated with an X-ray ring and torus. 
Objects known to display this morphology, likened by some authors to the shape of a cross- 
bow, include the Crab (Hester et al. 1995; Weisskopf et al. 2000; Bietenholz et al. 2001; 
Hester et al. 2002; Sollerman 2003), Vela (Pavlov et al. 2001), G320.4-1.2 (Gaensler et al. 
2002), G54.1±0.3 (Lu et al. 2002), 3C58 (Slane et al. 2002), SNR 0540-69.3 (Gotthelf & 
Wang 2000), and G0.9±0.1 (Gaensler et al. 2001), imaged variously by the Chandra X- 
Ray Observatory (CXO), Hubble Space Telescope (HST), Very Large Telescope (VLT), and 
Very Large Array (VLA), although there are counter-examples too, e.g. Gil. 2— 0.3 (Kaspi 
et al. 2001). The knots and wisps are interpreted, respectively, as the polar and equatorial 
zones of the magnetized, collisionless shock terminating the pulsar wind (Gallant & Arons 
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1994; Hester et al. 1995; Gaensler et al. 2002; Komissarov & Lyubarsky 2003; Spitkovsky 
& Arons 2004). The brightness asymmetries are ascribed to Doppler boosting. Ideas for 
explaining jet collimation and the observed jet-torus structure include magnetic hoop stress 
in the relativistic wind upstream from the shock (Begelman & Li 1994; Bogovalov 2001), the 
anisotropic energy flux in a force-free, monopole wind (Lyubarsky & Eichler 2001; Lyubarsky 
2002; Komissarov & Lyubarsky 2003; Del Zanna et al. 2004), magnetic hoop stress or hy- 
dromagnetic instabilities in the subsonic, downstream flow (Begelman & Li 1992; Begelman 
1998; Lyubarsky & Eichler 2001; Lyubarsky 2002; Melatos 2002, 2004), interactions with 
moving ejecta (Pavlov et al. 2003), and the helical structure of a wave-like, displacement- 
current-dominated wind (Usov 1994; Hester et al. 1995; Melatos & Melrose 1996; Melatos 
1998, 2002). 

The subarcsecond features of PWN are highly variable. In the Crab PWN, the optical 
wisps are seen to change brightness and position in time-lapse HST images every six days, 
receding from the pulsar concentrically in a wave-like pattern at 0.5c, while the optical 
knots jump around erratically on the same time-scale (Hester et al. 1995, 2002). In con- 
temporaneous CXO images, taken every 22 days, the X-ray ring inside the torus is resolved 
into more than 20 knots that brighten and fade irregularly, while the X-ray jet transports 
blobs of material and bow- wave-shaped structures radially outward at ~ 0.4c (Hester et al. 
2002; Mori et al. 2002). In time-lapse CXO images of the Vela PWN, Pavlov et al. (2001, 
2003) discovered changes of up to 30 per cent over several months in the brightnesses and 
spectra of the wisps, jet and centrifugal (0.3-0. 6c) knots, coherent bending of the outer jet 
over 16 days, and changes in knot brightness over just two days. The variability has been 
ascribed to an ion-cyclotron instability at the shock front (Spitkovsky & Arons 2004), to 
a synchrotron cooling instability (Hester et al. 2002), and to a nonlinear Kelvin-Helmholtz 
instability (Begelman 1999). 

In this paper, we report on the first near-infrared, adaptive-optics observations of the 
wisps and jet of the Crab PWN. The data offer high angular resolution (optimum 0."19, av- 
erage 0."23 in K'), high time resolution (0.6 ks, resolving the light-crossing and ion-cyclotron 
time-scales of the narrowest features for the first time), and the first feature-specific contin- 
uum color spectra extending from near-infrared to ultraviolet wavelengths. The observations 
are described in §2, the light curves and spectra of individual features are presented in §3 
and §4, and the results are interpreted physically in §5. 
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2. Observations 

The center of the Crab PWN was observed on 2002 February 6, 7, and 8 with the 
Hokupa'a/QUIRC adaptive optics instrument on the Gemini North Telescope (proposal 
GN-2002A-Q-16). Hokupa'a is a natural guide star, curvature sensing system with 36 ele- 
ments, coupled to a near-infrared camera, QUIRC, consisting of a 1024 x 1024 HgCdTe array 
with plate scale 0."0200, dark current < 0.1 electrons per second, and read-out noise 15-30 
electrons. The system is described in detail by Graves et al. (1998). Hokupa'a is capable 
of locking on to point sources with R < 15 mag nominally and R < 17 mag in practice. 
It is therefore ideal for imaging the Crab PWN, where there are two suitable guide stars 
within 4" of the center of the nebula (cf. isoplanatic radius ~ 30"): a field star (Star I) at 
a(2000) = 05 h 34 m 32.2 s and 5(2000) = ±22°00'57", with R = 14.8 ±0.2 mag, which we chose 
to use, and the Crab pulsar, with R = 16.3 mag pulse-averaged (Eikenberry et al. 1997). 
The wavelengths and bandpasses of the QUIRC J and K' filters are 1.25/zm/0.17/zm and 
2.12 /xm/0.41 jum respectively. We achieved resolutions of FWHM = 0."19-0."30 in K' and 
FWHM = 0. "24-0. "37 in J in this sequence of observations, where FWHM refers to the 
full- width half-maximum of the point spread function (PSF). The seeing fluctuated by up 
to ±0.05" in K' and ±0.08" in J over intervals of 0.6 ks. 

We obtained a sequence of 20 x 4 x 0.12 ks exposures spaced in a four-point, 3" dither 
pattern, which were subsequently combined into 20 frames, as follows: 5 x K' then 3 x J 
(2002 February 6), 3 x J then 5 x K' (2002 February 7), and 4 x K' (2002 February 8). 
Conditions were excellent on 2002 February 6. On the following nights, observations were 
affected by light cloud and wind. Data were gathered without interruption during all three 
nights, yielding a sampling time of 0.6 ks between frames (after read-out), with two brief 
exceptions: a five-point dither was accidentally performed on 2002 February 6, and two 
six-point dithers were required on 2002 February 8, when the telescope lost its guiding. All 
features were observed with good signal-to-noise. In a typical 0.12 ks exposure in K' , we 
accumulated w 4.9 x 10 2 counts/pixel for the pulsar and guide star, « 4.5 x 10 1 counts/pixel 
for the brightest extended feature (the sprite; see §3), and « 2.7 x 10 1 counts/pixel for the 
faintest extended feature (the faint wisp; see §3), after subtracting the nebula background 
(2.9-5.1 x 10 3 counts/pixel, or 1.1-1.9 Jyarcsec" 2 ). 

We assembled a data reduction pipeline in IRAF 1 to subtract bias, dark, and sky frames 
and divide by flat fields in the standard way. The three bright point sources in the field were 
excised from the sky frames with care, to avoid creating false shadows by oversubtraction. 
After trimming a border 20 pixels wide to remove faulty edge pixels, we median combined 
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each set of four dithered images, rejecting one high, then trimmed a border 3" wide to exclude 
the region where the four dithered exposures do not overlap. 

An important issue with any adaptive optics observation is the degree by which the 
PSF changes across the field of view. In this work, we can quantify the effect directly by 
examining the three point sources in the field (the guide star, the pulsar, and a field star 6" 
south of the pulsar, labeled Star II), which happen to be well separated. We find that the 
PSF is nearly axisymmetric at the southern field star, elongated in an east-west direction at 
the pulsar, and elongated in a northeast-southwest direction at the guide star, with minor 
and major axes in the ratio ~ 0.8 (although the isophotes are not strictly elliptical). This 
essentially precludes accurate photometry of the inner knot, located 0."6 from the pulsar, 
for the reasons set forth in §3.2, without prejudicing photometry of large, extended features 
like the wisps. Another issue is how the PSF varies as a function of time. We do not detect 
any change in the shape of the PSF at the locations of the three field stars when comparing 
isophotes from successive exposures by eye. However, there is indirect evidence that slight 
yet rapid changes do occur; we find that the total flux from the guide star within an aperture 
of radius 3.1 FWHM fluctuates by ±6 (±4) per cent in K' (J) over 0.6 ks after background 
subtraction (§3.1), accompanied by fluctuations of ±9 (±4) per cent in the southern field 
star. These changes feed into the measurement uncertainties calculated in §3 and §4, as the 
field stars serve as flux calibrators. 



3. Kilosecond variability 

Figure 1 shows the center of the Crab PWN at 0."19 resolution (FWHM) in K', as it 
appeared on 2002 February 6. Its crossbow- like morphology is evident. The subarcsecond 
features in the termination shock of the pulsar wind are cylindrically symmetric about the 
projected rotation axis and proper motion of the pulsar, determined from HST astrometry 
(Caraveo & Mignani 1999). 

1. The wisps are interpreted as shock structures in the equatorial plane of the pulsar 
wind (latitude ±10 deg), in the neighborhood of the X-ray ring and torus (Hester et al. 
1995; Weisskopf et al. 2000; Hester et al. 2002; Mori et al. 2002; Spitkovsky k Arons 
2004; Komissarov & Lyubarsky 2003; Del Zanna et al. 2004). The faint and bright 
wisps to the northwest of the pulsar, labeled in Figure 1, mark ion-driven magnetic 
compressions (at the first and second ion turning points) in the ion-cyclotron model 
of the shock (Gallant & Arons 1994; Spitkovsky & Arons 2004); they are probably 
analogous to the features labeled 5 and E by Gaensler et al. (2002) in another PWN, 
G320.4— 1.2. The position and brightness of these two wisps, several less prominent 
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wisps, and their fibrous substructure are known to change on a time-scale as short as 
six days in the optical (Hester et al. 2002) and 22 days in X-rays (Hester et al. 2002; 
Mori et al. 2002). 

2. The sprite can be interpreted as a polar shock, lying on the rotation axis (colatitude 
±15 deg) at or near the base of the polar X-ray jet (Weisskopf et al. 2000; Hester et al. 
2002; Mori et al. 2002). It can also be interpreted as a mid-latitude arch shock between 
the polar outflow and equatorial backflow in a pressure-confined, split-monopole neb- 
ula, situated at the tangent point of the line of sight (and therefore Doppler boosted) 
(Komissarov & Lyubarsky 2003; Del Zanna et al. 2004). Its shape, doughnut-like 
with a central rod (§3.2), changes irregularly on the same time-scale as the wisps; in 
CXO images, the sprite appears to be the launching point for blobs and 'bow waves' 
ejected along the X-ray jet (Hester et al. 2002; Mori et al. 2002), although these may 
also be unstable motions in the vicinity of the mid-latitude arch shock (Komissarov & 
Lyubarsky 2003). 

3. The inner knot is a barely resolved, flattened (§3.3) structure partly obscured in Figure 
1 by the pulsar PSF. It too can be interpreted as a polar feature, lying ps 5 times closer 
to the pulsar than the sprite. If the sprite marks the polar termination shock, the inner 
knot sits in the unshocked pulsar wind, and its physical origin is unknown (Hester et al. 
1995; Melatos 1998, 2002). If the sprite is a mid-latitude arch shock, the inner knot 
originates from a part of the arch shock nearer the base of the polar jet, which is 
pushed inward (relative to the wisps) because the energy flux in a split-monopole wind 
is lower at the poles than at the equator (Komissarov & Lyubarsky 2003). 

4. A conical halo is visible at intermediate latitudes, midway between the pulsar and faint 
wisp in Figure 1, the near-infrared counterpart of an optical feature noted in HST data 
by Hester et al. (1995). We do not discuss the halo further in this paper, as it is too 
faint for accurate near-infrared photometry. 

In this section, we examine the variability of four features — the bright wisp, faint wisp, 
sprite, and inner knot — in the near-infrared over time-scales as short as 0.6 ks, extending 
previous studies of the Crab PWN with HST (sampling time six days) and CXO (sampling 
time 22 days) (Hester et al. 2002; Mori et al. 2002). The first (and most challenging) step, 
subtracting the time-dependent nebula background, is discussed in §3.1. Light curves of the 
features are presented in §3.2 and §3.3. 
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3.1. Nebula background 

It is difficult to characterize and hence subtract the background in Figure 1, because 
there is no unique way to disentangle the contributions from the nebula and sky, given 
that the nebula is ubiquitous, nonuniform, and time-dependent. The surface brightness 7 G 
observed in any pixel is the sum of flux from the nebula background (7 n ) and any feature 
(If) occupying that pixel, with 7 G = (7 n + I{)i] a + 7 a , where r\ a and 7 a denote the atmospheric 
absorption coefficient and sky brightness respectively. There are two problems in extracting 
If from 7 G . First, there are no pixels empty of both features and nebula, so we cannot 
measure 7 a directly. In our analysis, we assume that r\ a and l a are uniform across the field of 
view, but both parameters vary markedly from exposure to exposure, as quantified below. 
Second, 7 n fluctuates stochastically from pixel to pixel, so we cannot estimate 7 n behind a 
feature by interpolating 7 n directly from neighboring, feature-free pixels. 

To overcome these problems, we average the observed brightness of all feature-free pixels 
(I' ) in the field to obtain I n i] a + I a — (I' Q ) and hence, approximately, Ifrj a — 7 — (I' Q ) along 
any line of sight with a feature, under the assumption (justified below) that there is no large- 
scale gradient of nebula brightness across the field. The uncertainty in this estimate of Ifi] a is 
given by the width of the I' Q distribution, measured below. Without extra, exposure-specific 
information, it is impossible to determine r\ a and 7 f independently. However, we are interested 
here in the variability of features rather than their absolute brightness. Consequently, we 
can normalize the observed brightness of any feature, J ,f, to the observed flux 7 0jg of an 
intrinsically steady point source J g (e.g. the guide star or pulsar), after subtracting the 
nebula background, to obtain 7f/7 g = (I Q ,i — (I' ))/(Io,g ~ (-0)- If ^ were necessary to 
determine If absolutely, we would need to measure 7 g independently, e.g. from J and K' 
photometry of the the Crab pulsar (Eikenberry et al. 1997). 

How accurate is the above approach? A histogram of pixel counts for a single K' frame, 
excluding point sources and extended features, is presented in Figure 2a (solid curve). The 
mean, median, and standard deviation (51) of the distribution in Figure 2a are 5033, 5024, 
and 13.2 counts respectively. The distribution is not Gaussian, cutting off sharply at ±3.5(57, 
and it is narrower than Poisson (51 = 0.28(7q) 1 / 2 ) because 7 n is correlated in neighboring 
pixels. For our 20 images, we find 51 in the range 10-30 counts, while the median varies 
markedly in the range 2.9-5.1 kcounts. These statistics are corroborated by the dashed 
curve in Figure 2a, a histogram of pixel counts averaged over 50 x 50 pixel blocks (chosen 
to roughly match the dimensions of the knot-like features of interest). By inspecting the 
image directly, we isolate the blocks that appear to be empty of features, obtaining a mean, 
median and standard deviation of 5033, 5024, and 13.4 counts per block, in close agreement 
with the solid curve (they are indistinguishable to the eye). We also verify by inspection that 
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there is no large-scale gradient in counts per block across the field of view, confirming that 
r? a and J a are uniform within the statistical uncertainty 51. Finally, in Figure 2b, we present 
the histogram of pixel counts in an annular aperture of radius 6.0 FWHM centered on the 
guide star. Annular and polygonal apertures are provided for background measurements in 
the Gemini IRAF software and are employed in §3.2 and §3.3. The statistics are consistent 
with Figure 2a. We find that the median background in the annulus differs by at most 8 
counts from the median of the field in all 20 images, well within the standard deviation 51, 
and is arguably a more accurate estimate of the background locally. 

We estimate the uncertainty in our fluxes as follows. The absolute uncertainty in hr] a , 
the total flux minus the background, is given by (e\ + e 2 , + e 2 ) 1 / 2 , where e\ is the square 
root of the counts after background subtraction (corrected for the ADU-photon ratio), e 2 
equals N^ 2 5I (N { is the number of pixels in the aperture enclosing the feature), and e 3 
is given by (Nf/N s ) 1 / 2 e2 (N s is the number of pixels in the aperture estimating the sky). 
A similar uncertainty attaches to the guide star I g i] a . Note that e\ represents the Poisson 
fluctuation in the intrinsic flux of the feature; e 2 measures the uncertainty in the background 
contribution to the total flux, characterized by Figure 2a and 51 (not Poissonian); and 
is the uncertainty in the background level subtracted from the total flux, corrected for the 
relative sizes of the sky and feature apertures. 



3.2. Termination shock: wisps and sprite 

In this section, we examine the variability of the equatorial and polar zones (wisps and 
sprite) of the termination shock in the Crab PWN. 

Figure 3 shows two enlarged images of the bright and faint wisps in J band, taken 1.2 ks 
apart on 2002 February 7, after nebula subtraction and normalization to the guide star (Ii/I g ; 
see §3.1). The frames should be identical if there is no change in the intrinsic brightness 
of the wisps, yet differences between them are readily apparent (although the features are 
not displaced). To quantify the changes, we use the polymark tool in IRAF to specify a 
polygonal aperture enclosing each wisp, as drawn in Figure 3. The fluxes enclosed by the 
apertures, after nebula subtraction and normalization, are plotted as functions of time in 
Figure 4. J- and i^'-band data are both displayed; uncertainties are calculated according to 
the recipe in §3.1. For the bright wisp, we find maximum flux changes (relative to the mean 
level) of ±24 ± 4 per cent in J and ±12 ± 7 per cent in K', occurring in the space of 1.2 ks 
on 2002 February 7, and smaller fluctuations at other times. Moreover, the light curves of 
the bright and faint wisps appear correlated to some degree in both filters. The detection 
of variability is marginal in K' but more statistically significant in J, where the nebula 
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background is less prominent. We find, by experimentation, that the results are essentially 
independent of the choice of aperture, while the time- dependent PSF has a minimal effect 
on the measured flux of the extended features (for the guide star, the effect is included in 
the measurement uncertainty; see §2). Nevertheless, new observations — preferably by an 
independent party using a different instrument — need to be made before variability of the 
wisps on such short time-scales can confidently be claimed. 

Figure 5 shows two enlarged images of the sprite in J band, taken on 2002 February 
6 and 8. It is interesting to note its doughnut-like structure, symmetric about the pulsar's 
rotation axis, as well as the short, bent rod emerging from its center, seen clearly here for 
the first time and corroborating the observation by Hester et al. (2002) that the sprite is 
often center-filled (especially in X-rays). Knot-like features in another PWN, G320.4— 1.2, 
numbered 2 and 3 by Gaensler et al. (2002), may be analogs of the sprite. 

In common with the wisps, there are visible differences in the brightness (but not the 
position) of the sprite in the two images, even after nebula subtraction and normalization - 
not just the absolute brightness, but also, more significantly, the brightness contrast between 
the rod and doughnut, which is less likely to be affected by PSF nonuniformity and imperfect 
sky/nebula subtraction. (The flux of the guide star is equal to within 4 per cent in the two 
snapshots.) To quantify the brightness changes, we define apertures enclosing the rod and 
the whole sprite, and plot the nebula-subtracted, normalized aperture flux If/I g versus time 
in Figure 6. We measure the following maximum flux changes (relative to the mean): ±13±4 
per cent (rod, J), ±15 ± 7 per cent (rod, K'), ±9 ±4 per cent (whole sprite, J), and ±8 ± 7 
per cent (whole sprite, K'). As for the wisps, the detection of variability is marginal in K' 
and somewhat more significant in J, especially for the rod, as seen in Figure 5. 



3.3. Pulsar wind: inner knot 

The inner knot, discovered by Hester et al. (1995) in optical HST data, is displaced 
0."65 from the pulsar along the axis of symmetry of the PWN, and is resolved by HST to be 
m 0."2 thick (Hester et al. 1995). No counterpart has been detected unambiguously at X-ray 
wavelengths, although there is a hint of a southeasterly 'bump' protruding from the pulsar 
in CXO images, e.g. in Figure 5 of Hester et al. (2002). It is also possible that an analogous 
feature, named feature 1, has been discovered in CXO images of another PWN, G320.4— 1.2 
(Gaensler et al. 2002). In Figure 7, we present a brightness map of the inner knot, after 
subtraction of the PSF. It is clear, from the isophotes (solid contours) in particular, that the 
feature is flattened, not spherical, although it is hard to discern its shape exactly because it 
is barely resolved in our highest-resolution K' data and the PSF subtraction is imperfect. 
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Photometry of the inner knot is complicated by its proximity to the pulsar, whose flux 
contaminates the knot unpredictably from image to image as the seeing fluctuates. The 
standard approach, modeling and subtracting the pulsar PSF, is attempted in Figure 7, 
but the result is unreliable; see §2. Faced with these difficulties, we test for variability 
of the inner knot without measuring its brightness directly in two complementary ways. 
In the first test, we measure the fluxes ip(ri) and Ipfa) (after subtracting (I' )) enclosed 
by two circular apertures centered on the pulsar, of radii r x = 0.75 FWHM (including as 
much PSF as possible but excluding most of the knot) and r 2 = 5.0 FWHM (including 
the PSF and knot). The apertures are drawn in Figure 7. The ratio of these fluxes, after 
nebula subtraction, would be the same as the ratio of the fluxes I s (ri) and I g (r 2 ) enclosed by 
identical apertures around the guide star if there were no inner knot, because the cylindrically 
averaged PSF is uniform to within ±6 (±4) per cent in K' (J) across the field of view (see 
§2). Therefore, the flux difference Ik/ I g = 7 p (r 2 )// g (r 2 ) — Ip{ r i) / Ig{ r i) can be attributed 
to the presence of the inner knot, and any change in Ik/Ig from image to image is evidence 
that the inner knot varies intrinsically. In Figure 8, we plot Ik/Ig as a function of time for 
the full sequence of observations. The data are consistent with no variability, within the 
measurement uncertainties. For example, in the K' band, we find maximum peak-to-peak 
changes in J k /J g of 0.045 ± 0.060 over 1.2 ks and 0.054 ± 0.063 over 48 hours. 

A second test provides a consistency check: we measure directly the flux enclosed by a 
circular aperture of radius 10 pixels, centered on the knot. The flux thus measured includes 
leakage from the pulsar PSF, the amount of which varies from image to image along with the 
seeing, as noted above, but strong intrinsic variations in knot brightness could still overwhelm 
this effect. In Figure 9, we plot the flux in the 10-pixel aperture as a function of time, after 
nebula subtraction and guide star normalization. The result is consistent with Figure 8: 
there is no significant detection of variability within the measurement uncertainties, with a 
maximum peak-to-peak change of 0.011 ± 0.007 over 1.2 ks on 2002 February 7. 

On the strength of the data presented here, we are unable to say whether or not the 
inner knot is variable on time-scales of 0.6 ks to 48 hours. New observations — preferably 
by an independent party using a different instrument — are required to settle the issue, and 
dedicated PSF calibration frames will be essential if adaptive optics are used. 
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4. Feature-specific color spectra 

4.1. Near infrared: wisps, sprite, and inner knot 

In this section, we measure the J-to-K' color spectra of the faint wisp, bright wisp, 
sprite, rod, and inner knot. As these features may vary on time-scales shorter than the 
minimum interval between exposures, we sum the 14 K' images and 6 J images in our data 
set to obtain time-averaged spectra. 

To the best of our knowledge, calibrated J-to-K' spectra of Stars I and II have not 
been published. Therefore, to calibrate the spectra of the various subarcsecond features in 
the PWN, we are forced to redo the photometry in §3 by normalizing nebula-subtracted 
fluxes with respect to the pulsar, whose phase-averaged, near-infrared, color spectrum was 
determined by Eikenberry et al. (1997) using the Solid State Photomultiplier on the Multiple 
Mirror Telescope. After dereddening, the pulsar's spectrum is fairly flat in the near infrared 
(F u oc z/ ' 36 ); see the third row of Table 1. Note that the fluxes in Eikenberry et al. (1997) 
include the inner knot (unresolved from the ground). By comparing with the total flux of the 
pulsar plus knot in our images, we derive a zero magnitude reference for calibrated aperture 
photometry of the extended features (e.g. wisps, sprite, and rod). 

To measure the flux of the knot, we remove the contribution from the pulsar by scaling 
it to the azimuthally averaged PSFs of Stars I and II, such that the normalization is the 
same at a radius of 10 pixels. The combined profile of the pulsar and knot is found to match 
the stellar PSFs up to a radius of w 20 pixels (to within 0.05 mag in K' and 0.03 mag in J) 
but deviates beyond due to the knot excess. We measure the flux difference between radii 
of 20 and 60 pixels (where the profile merges into the background) and adjust the results in 
Eikenberry et al. (1997) to give the calibrated knot flux. Note that the differential method 
of detecting the knot (1^/ I g ) employed in §3.3 does not yield a calibrated flux. 

Time-averaged J- and i^'-band fluxes are presented in Table 1 for the point-like and 
extended subarcsecond features identified in Figure 1, together with the spectral index a 
of each feature assuming its flux density scales as F v oc v a (as for synchrotron radiation, 
although not necessarily for all nonthermal processes, e.g. synchro- Compton radiation; see 
§5). The fluxes for the extended features are quoted per unit length, although it is possible 
that we are marginally resolving the wisps across their width. The field stars are included for 
reference, as their spectra have not been published previously. Three sources of uncertainty, 
summed in quadrature, contribute in Table 1. First, the central wavelength of the K filter 
used by Eikenberry et al. (1997) is 0.08 /zm greater than for the QUIRC K' filter. Second, 
there is scatter in the 20-pixel flux and 20-to-60-pixel offset. Third, the extinction corrections 
are uncertain (0.63 ± 0.03 mag in J and 0.81 ± 0.03 mag in K'). 
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Several interesting conclusions emerge from these data. First, the inner knot is clearly 
the reddest feature in the region, with a ~ —0.8. This is understandable if the inner knot is 
produced in the pulsar wind by a different radiation mechanism than other features. However, 
it is surprising if the inner knot is part of the same arch shock that produces the sprite. 
Second, the polar or mid-latitude sprite and rod have flatter spectra (a = —0.21 ± 0.13) than 
the equatorial wisps (a = —0.56 ±0.12). Yet all these features are synchrotron emitting 
elements of the same termination shock, albeit at different latitudes. Third, there is no 
feature in the region whose spectrum matches smoothly from the near infrared to X-rays. 
The theoretical implications of these results are considered further in §5. 



4.2. Ultraviolet: inner knot 

We extend the spectrum of the inner knot in Table 1 with an independent measurement 
of the ultraviolet flux of this feature. Gull et al. (1998) observed the Crab pulsar with the 
HST Space Telescope Imaging Spectrograph (STIS) NUV-MAMA detector on 1997 August 
7 through the low dispersion G230L grating. The observations were made using a 2" x 2" 
aperture which included the inner knot, with part of the exposure (2 x 2.4 ks) in TIME- 
TAG mode. We acquired the archival data, barycentered the photons using standard STIS 
routines, and extracted a 'slit' image using photons from 0.25 rad of phase spanning the 
pulse minimum, thereby gating out the pulsar. Approximately two per cent of the unpulsed 
flux remained, but, after subtracting a scaled version of the on-pulse PSF, a clear excess was 
found at the projected position of the inner knot 0.3 ±0.2" on one side of the pulsar. The 
intensity profile agrees with that expected from a one-dimensional collapse of direct HST 
images. Assuming a flat a as spectral index over the NUV band (160-320 nm), we find a 
summed inner knot flux (3.3±0.1) x 10~ 2 times the unpulsed flux of the pulsar. De-reddening 
with the best fit value E(B - V) = 0.52 (Sollerman et al. 2000) yields F v = 0.12 ± 0.01 mJy 
for wavelengths in the range 0.16-0.32 /xm. Note that the ultraviolet flux is consistent with 
the near-infrared spectrum measured in §4.1, which extrapolates to give F u = 0.11 mJy at 
0.32 /im for a = -0.8. 



5. Discussion 

In this paper, we report on the first near-infrared, adaptive-optics observations of the 
wisps and jet of the Crab PWN, comprising 20 J- and fT'-band snapshots taken at 0."19- 
0."37 resolution with the Hokupa'a/QUIRC camera on the Gemini North Telescope. The 
data contain tantalizing — albeit inconclusive — evidence that subarcsecond features in 
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the termination shock of the Crab PWN vary intrinsically in J-band brightness by ±24 ± 4 
(wisps) and ±14 ±4 (sprite) per cent on time-scales as short as 1.2 ks. The principal sources 
of uncertainty are the nonuniform, unsteady nebula background and PSF. The data also 
suggest that the near-infrared spectra of polar features in the termination shock are flatter 
(e.g. sprite, F v oc //-°- 21±0 - 13 ) than the spectra of the equatorial wisps (F u oc //-°- 56±0 - 12 ) ) 
except for the steep-spectrum inner knot (F v oc z/" 08 ), which may lie in the unshocked 
pulsar wind. This result is supported by an independent measurement of the ultraviolet flux 
of the inner knot, obtained by reanalyzing archival, time-tagged, HST STIS data. 



5.1. Ion cyclotron and fire hose instabilities 

Why, physically, might the nebula vary on time-scales as short as 1.2 ks? According 
to one hypothesis, modeled numerically by Spitkovsky & Arons (2004), the wind contains 
ions which drive a relativistic cyclotron instability at the termination shock. The instability 
exhibits limit cycle dynamics: ion bunches and magnetic compressions are launched down- 
stream periodically at roughly half the ion-cyclotron period, |Tj = nAiUip'-yi/ZieB, where A4 
and Zi are the atomic number and charge, 7» is the preshock Lorentz factor, and B is the 
postshock magnetic field. In order to fit the separation of the innermost wisps, one must take 
T; ps 3.2 x 10 7 s, consistent with 7; = 7.8 x 10 5 and B — 16 fiG at a radial distance r = 0.1 pc 
from the pulsar (for A^jZ\ = 1) (Spitkovsky & Arons 2004). Faster variability is expected 
nearer the pulsar, because the magnetic field in the wind scales as B = 16(r/0.1 pc) -1 /iG. 2 
However, we find 7] ps 2 x 10 7 s and 7] ~ 2 x 10 6 s at the sprite (r ps 0.05 pc) and inner knot 
(r ps 0.007 pc) respectively, slower than the variability observed. Relativistic Doppler boost- 
ing does not improve the agreement; the time-scale oc 7f 2 is unchanged downstream (7; ~ 1) 
and too short upstream (7^ ~ 10 6 ). The electron-cyclotron period, T e = {m e Zi/ m v Ai)T u does 
fall in the observed range, but it is hard to see how to maintain coherent limit-cycle dynamics 
in the electrons when they are randomized rapidly at the shock by ion-driven magnetosonic 
waves (Gallant & Arons 1994; Spitkovsky & Arons 2004). We are therefore inclined to rule 
out a cyclotron origin of the observed kilosecond variability in the near infrared. 

The argument against a cyclotron origin of the kilosecond variability assumes that all 
the energy in the unstable (compressional) ion-cyclotron-magnetosonic waves resides in the 
fundamental. This need not be so. The frequency spectrum of the waves is quite flat in 
one-dimensional simulations (Hoshino et al. 1992); significant power is deposited at high 



2 Fastcr variability is also expected at certain special phases in the month-long ion cycle, e.g. in the 
neighborhood of a moving wisp, where the plasma is stirred up (A. Spitkovsky, private communication). 
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harmonics (up to orders ~ m p /m c ) provided that parametric three-wave decays do not 
destroy the coherence of the waves, a plausible concern in a realistic, three-dimensional 
plasma (J. Arons, private communication). 

Another possible mechanism, applicable especially to the knots in the polar jet, is the 
relativistic fire hose instability, driven by anisotropy of the kinetic pressure parallel (Py) 
and perpendicular (Pj_) to the magnetic field B. Noerdlinger & Yui (1969) showed that the 
growth time in the bulk frame of the jet (denoted by primes) is given by 

1 _ eg 1.3-i?' 2 /M^-n)] 

Tfh m e < 7 ') (l + cVO^ + S ' U 

where (7') is the thermal Lorentz factor, v® = B' 2 /p p' is the Alfven speed, and p' is the 
density of the cold ion background. As long as the condition Py — P± > 0.77B' 2 /p is met, 
the minimum growth time is given by T{ hmin = 0.14m e (7 / )/ei? / for v' A ^> c. Upon Lorentz 
transforming to the observer's frame, we obtain (i) Tfh jm i n = 0.14m c (7)/eP if B is radial 
in the jet, or (ii) Tfh, m m = 0.147;m c (7)/eP if B is helical or toroidal in the jet, assuming 
(7} m 71(77. (The validity of the last assumption depends subtly on the exact form of 
the electron distribution.) Immediately upstream from the termination shock of the jet, we 
have (7) = 1, B = 2.7 pG, and 7 = 7.8 x 10 5 , implying T fllimin = 2.3 ks in scenario (ii). 
Immediately downstream from the termination shock, we have (7) « 1 x 10 6 , B = 8 pG, 
and 71 = 1.1, implying Tfh, m in ~ 1-1 ks in scenarios (i) and (ii) — intriguingly close to the 
observed time-scale. Note that the fire hose instability only occurs for P|j > P' ± ; in the 
reverse situation, a mirror instability exists for P' ± — P|| > 8.9B' 2 /p , with growth time 
~ Tqj. Gallant & Arons (1994) argued that P!JP'j_ increases from zero to unity downstream 
from the shock — the adiabatic index (3 + P!JP' ± )/(2 + Py/Pj.) decreases from 3/2 to 4/3 
as pitch-angle scattering isotropizes the electrons — but Pj > P' ± cannot be excluded. 

Despite appearances, it is unlikely that the fire hose instability causes the unsteady, 
serpentine motions observed in the Vela X-ray jet on time-scales between one day and several 
weeks (Pavlov et al. 2003), because the growth time appears to be too short. We estimate 
Pfh,min = 0.147im c (7)/eP « 0.1((j/10" 3 ) 1/2 ks, taking (7) = 7 x 10 4 [pair multiplicity w 10 3 ; 
see Figure 17 of Hibschman & Arons (2001)] and B = 0.15a 1 / 2 mG (radial or toroidal), where 
a is the ratio of Poynting to kinetic energy flux at the base of the jet. Pavlov et al. (2003) 
suggests an alternative scenario, in which the end of the jet is bent by an external wind 
while the knots in the body of the jet are produced by hydromagnetic (kink and sausage) 
instabilities on the local Alfven time-scale. 

Our time-lapse observations resolve the light-crossing time-scale of the smallest features 
in the field, e.g. 2.2 days for 0."19 at 2.0 kpc. Therefore, if the kilosecond variability we 
observe is real, it must arise from (i) a pattern traveling at a superluminal phase speed, or 
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(ii) relativistic Doppler boosting in the upstream collimated outflow (7 ; ~ 2 times the light- 
crossing time-scale; cf. millisecond variability of unresolved gamma-ray bursters). Our results 
are consistent with previous observations that also detected significant variations on, or faster 
than, the light-crossing time-scale: Hester et al. (2002) observed the optical knots in the Crab 
PWN to vary over six days, and Pavlov et al. (2003) observed the X-ray jet in the Vela PWN 
to vary over just two days. 



5.2. Radiation mechanisms 



The near-infrared spectral indices displayed in Table 1 are curious in several respects. 
First, the inner knot has a steeper spectrum than every other feature in the region — not just 
in our data, where the J-band flux is uncertain to ±50 per cent, but also in data obtained 
with the Infrared Spectrometer And Array Camera on the VLT in 0. "65-0. "88 natural seeing 
(Sollerman 2003). One explanation is that the inner knot is physically different from the 
other features: the wisps and sprite are part of the termination shock and emit synchrotron 
radiation, e.g. from ion-cyclotron- heated electrons (Gallant & Arons 1994), whereas the inner 
knot lies upstream in the unshocked pulsar wind and emits synchro-Compton radiation, e.g. 
from electrons heated by magnetic reconnection (Coroniti 1990; Lyubarsky & Kirk 2001) or 
parametric instabilities (Melatos & Melrose 1996; Melatos 1998, 2002) in a wave-like wind. 
This explanation conflicts with recent simulations which suggest that the inner knot is syn- 
chrotron emission from an arch shock between the polar outflow and equatorial backflow in 
the nebula (Komissarov & Lyubarsky 2003). Unfortunately, we cannot discriminate between 
these two possibilities spectrally, because both synchrotron and synchro-Compton radiation 
yield a — (1 — p)/2 given a power-law electron distribution N("f) oc 7~ p (Blandford 1972; 
Leubner 1982). [Monoenergetic electrons in a large-amplitude wave emit an inverse Comp- 
ton spectrum F v oc v at frequencies below 0.l7 2 a 3 f2 ~ 10 16 Hz, where a = 1-10 is the wave 
nonlinearity parameter (Melatos & Melrose 1996; Melatos 1998) and Q is the pulsar spin 
frequency, but this low-frequency tail is modified to F v oc z/( 1_p )/ 2 for power-law electrons.] 
Instead, we propose that the near-infrared polarizations of the inner knot and the sprite be 
measured. Theory predicts that synchrotron and synchro-Compton radiation are polarized 
perpendicular and parallel to the magnetic field respectively; one expects 50-80 per cent 
linear polarization from a linearly polarized large-amplitude wave if 0.5 < p < 5 (Blandford 
1972). Therefore, assuming that the magnetic geometry of the wave-like wind is similar at 
the inner knot and sprite (e.g. helical at high latitude), the polarization vectors of the two 
features are expected to be perpendicular if the inner knot originates from the unshocked 
pulsar wind and parallel if it is an arch shock. Moreover, if the inner knot originates from 
the unshocked wind, it should also be circularly polarized (degree ~ a -1 , independent of v). 
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A second curious property of the spectra in Table 1 is that the bright and faint wisps are 
steeper than the sprite and rod. In this case, there is no doubt that both sets of features are 
shock-related synchrotron emission at equatorial and polar latitudes respectively, whether 
the sprite is located at the working surface of the polar jet (Hester et al. 2002) or along 
the arch shock (Komissarov & Lyubarsky 2003; Del Zanna et al. 2004). So why do the 
spectra differ? One possibility is that the polarization of the large-amplitude wave in the 
wind zone, which is linear at the equator and circular at the pole, affects the acceleration 
physics in the shock ponderomotively; tentative indications to this end are emerging from 
recent particle-in-cell simulations (O. Skjaeraasen, private communication). 

Table 1 raises a third puzzle: the near-infrared spectra of some features do not extrap- 
olate smoothly to optical and X-ray wavelengths, while others do. For example, the HST 
V-band surface brightness of the bright wisp and the V^-band flux of the sprite are measured 
to be 49/xJy arcsec -1 and 4.1 /iJy respectively (Hester et al. 1995); the same quantities, ex- 
trapolated from Table 1, are predicted to be ~ 42/zJy arcsec -1 and w 15 yuJy respectively. 
VLT and HST observations show that the spectrum of the inner knot extends smoothly from 
near-infrared to optical wavelengths, with a = —0.75 ± 0.15 (Sollerman 2003), and this is 
corroborated (within larger uncertainties) by our data. On the other hand, extrapolation of 
the VLT spectrum of the inner knot to 0.16-0.32 /im yields F v = 0.35±0.05 mJy, significantly 
greater than the ultraviolet flux measured in §4.2. 

The spectrum of the bright and faint wisps is significantly steeper in X-rays than in 
the near infrared, with «x ~ —1-3 (Weisskopf et al. 2000) and q;nir — «x ~ 0.7, compa- 
rable to the steepening expected from synchrotron cooling (Y. Lyubarsky, private com- 
munication). The synchrotron cooling time for near-infrared-emitting electrons, t coo \ = 
2.8 x 10 4 (S/16^G)- 3 / 2 (i//2.4 x 10 14 Hz)" 1 ^ yr , is much longer than the flow time t flow 
across the wisps. On the other hand, the near-infrared spectral index of the sprite and 
rod (a = —0.21 ± 0.13) is significantly shallower than the wisps and similar to the average 
radio spectral index of the nebula (Y. Lyubarsky, private communication). This is either 
coincidental or highly surprising. The sprite and rod are shock features in which one has 
^flow *C i C ooi f° r near-infrared-emitting electrons and kilosecond variability is observed; they 
reflect electron acceleration at the present time. The radio electrons reflect the history of 
electron acceleration over the lifetime of the nebula and should be unaffected by the present 
dynamics of the sprite. This paradox is encountered in a related context: Bietenholz et al. 
(2001) observed that the radio and optical wisps travel radially in concert and display coor- 
dinated spectral index variations. 

We conclude by reiterating that the measurement uncertainties in our observations are 
substantial, due to the nonuniformity of the nebula and PSF. Consequently, the evidence for 
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kilosecond variability and spectral differences, while tantalizing, is inconclusive. Improved 
observations - - preferably by an independent party using a different instrument - - are 
essential to clarify the situation. 
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Fig. 1. — Near- infrared image of the central 20" x 20" region of the Crab PWN, constructed 
from four dithered, 0.12 ks exposures in the K' filter, taken on 2002 February 6 with a 
resolution of FWHM = 0."19. The principal features of the pulsar wind termination shock 
and its environment are labeled, revealing the 'crossbow' morphology characteristic of many 
PWN: arc-like, equatorial wisps and jet-like, polar knots symmetric about the projected 
rotation axis of the pulsar. Star I is the Hokupa'a guide star. Star II is another field star. 



4900 5000 5100 5200 4900 5000 5100 5200 

counts counts 

Fig. 2. — Statistics of the nebula background, a. Left. Frequency histogram of pixel counts 
in a typical K' image constructed from 4 x 0.12ks dithered exposures on 2002 February 
6, excluding point sources and extended features. Raw counts (solid curve) and counts 
averaged over 50 x 50 pixel blocks (dashed curve) are shown; the curves overlap and are 
indistinguishable to the eye. b. Right. Frequency histogram of pixel counts in an annulus, 
radius 6.0FWHM and one pixel thick, centered on the guide star (Star I). 
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Fig. 3. — Enlarged J-band snapshots of the bright and faint wisps, taken on 2002 February 
7, separated by 1.2 ks. Orientation is as for Figure 1, with north up and east left. Each image 
has been background-subtracted and normalized to the flux of Star I. The gray scale is the 
same in both images and is indicated by the central bar. The units on the bar indicate the 
normalized flux (defined as in Figure 4) in a 1" x 1" box. The apertures used to measure the 
total flux of the wisps, plotted in Figure 4, are drawn as dotted curves. The seeing and the 
guide star flux are the same in the two exposures to within 0."04 and 4 per cent respectively. 
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Fig. 4. — Light curves of the bright and faint wisps, after nebula subtraction and guide 
star normalization, as functions of time (in seconds) during the three nights of observations 



Fig. 5. — Enlarged J-band snapshots of the sprite (the fainter, larger feature on the left) and 
the rod (the brighter, smaller feature at center right), taken on 2002 February 7, separated 
by 1.2 ks. Orientation is as for Figure 1, with north up and east left. Each image has been 
background-subtracted and normalized to the flux of Star I. The gray scale is the same in 
both images and is indicated by the central bar. The units on the bar indicate the normalized 
flux (defined as in Figure 4) in a 1" x 1" box. The apertures used to measure the total fluxes 
of the sprite plus rod (plotted in Figure 6) and the rod alone (quoted in Table 1) are drawn 
as dotted polygons. The seeing and the guide star flux are the same in the two exposures to 
within 0."04 and 4 per cent respectively. 
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Fig. 6. — Light curve of the sprite, after nebula subtraction and guide star normalization, 
as a function of time (in seconds) during the three nights of observations (demarcated by 
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Fig. 7. — Enlarged K' image of the inner knot after subtraction of the pulsar PSF, plotted 
as a gray scale image with isophotes (solid contours) overlaid. The subtraction is imperfect 
because the PSF template is taken from the guide star, yet the PSF varies across the field 
of view (§2). Three apertures used in our photometry of the inner knot are also shown: two 
circles of radius 0.75 FWHM and 5.0 FWHM, centered on the pulsar, and a circle of radius 
10 pixels, centered on the brightness centroid of the inner knot. An annular region of inner 
radius 6.0 FWHM and one pixel width, also centered on the inner knot and used to estimate 
the background (§3.3), lies outside the borders of the image. 
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Fig. 8. — Light curve of the inner knot, plotted as the flux ratio difference J k / I g = 
/p(r 2 )//g(r 2 ) - Jp(ri)// g (r 1 ) for apertures with n = 0.75 FWHM and r 2 = 5.0 FWHM. 
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Fig. 9. — Light curve of the inner knot, plotted as the flux enclosed by a circular aperture 
of radius 10 pixels (after nebula subtraction and guide star normalization). This aperture 
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Point source 


J flux 
[mJy] 


K' flux 
[mJy] 


a 


length 

[arcsec] 


Star I 


6.25 (0.12) 


3.09 (0.07) 


1.24 




Star II 


1.01 (0.02) 


0.76 (0.02) 


0.50 




Pulsar 


3.21 (0.15) 


2.62 (0.17) 


0.36 




Knot 


0.32 (0.16) 


0.49 (0.17) 


-0.78 




Extended feature 


J flux 


K' flux 


a 


length 




[/xJy arcsec -1 ] 


[/xJy arcsec -1 ] 




[arcsec] 


Bright wisp 


69.4 (1.4) 


97.7 (2.2) 


-0.60 


7.3 


Faint wisp 


59.8 (1.3) 


80.3 (1.8) 


-0.52 


3.3 


Rod 


62.8 (1.3) 


68.7 (1.6) 


-0.16 


1.2 


Sprite 


46.2 (1.0) 


53.5 (1.2) 


-0.26 


2.5 



Table 1: De-reddened fluxes and power-law spectral indices (F„ oc v a ) for features present in 
the images, collected into point sources (including the inner knot; top) and extended features 
(bottom). Fluxes for the extended features are quoted per arcsecond of length; the lengths 
are shown in the right-hand column. Errors (la) are quoted in parentheses after the fluxes. 
A correction has been made for interstellar extinction (Eikenberry et al. 1997). 



